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ABSTRACT 

We present the number counts of Ks -band selected high redshift galaxy popu- 
lations such as extremely red objects (EROs), B-, z- & K -band selected galaxies 
(BzKs) and distant red galaxies (DRGs) in the AKARI NEP field. These high 
redshift galaxy samples are extracted from a multicolor catalog combining optical 
data from Suprime-Cam on the 8.2 m Subaru telescope with near-infrared data 
from the Florida Multi-object Imaging Near-IR Grism Observational Spectrom- 
eter on the Kitt Peak National Observatory 2.1 m telescope over 540 arcmin 2 in 
the NEP region field. The final catalogue contains 308 EROs (K s < 19.0; 54% 
are dusty star- forming EROs and the rest are passive old EROs), 137 star-forming 
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BzKs and 38 passive old BzKs (K s < 19.0) and 64 DRGs (K s < 18.6). We also 
produce individual component source counts for both the dusty star-forming and 
passive populations. 

We compare the observed number counts of the high redshift passively evolv- 
ing galaxy population with a backward pure luminosity evolution (PLE) model 
allowing different degrees of number density evolution. We find that the PLE 
model without density evolution fails to explain the observed counts at faint mag- 
nitudes, while the model incorporating negative density evolution is consistent 
with the observed counts of the passively evolving population. 

We also compare our observed counts of dusty star-forming EROs with a phe- 
nomenological evolutionary model postulating that the near-infrared EROs can 
be explained by the source densities of the far-infrared - submillimetre popula- 
tions. Our model predicts that the dusty ERO source counts can be explained 
assuming a 25 percent contribution of submillimetre star-forming galaxies with 
the majority of brighter Ks -band detected dusty EROs having luminous (rather 
than HR10 type ultra-luminous) submillimetre counterparts. We propose that 
the fainter Ks >19.5 population is dominated by the sub-millijansky submil- 
limetre population. We also predict a turnover in in dusty ERO counts around 
19< K s <20. 

Subject headings: galaxies: evolution — galaxies: formation — galaxies: fun- 
damental parameters (classification, colors, number counts) — galaxies: high- 
redshift — galaxies: starburst — galaxies: statistics — infrared: galaxies 



Introduction 



A quantification of the number density evolution of each galaxy population in successive 
epochs would provide stringent constraints on galaxy formation models. This decade has 
seen enormous progress in this field, particularly in the K -band selected high redshift 



galax y populations such as the extremely red objects (EROs, e.g. lElston. Rieke fc Rieke 



19881 . sometimes also referred to a s extremely red galaxies; ERG s), B-, z- and K -band 



selected galaxies (IDaddi et al.ll2004l . BzKs) and distant red galaxies (Ivan Dokkum et al.ll2003 
Franx et alil2003l . DRGsl 



EROs are selected to have very red optical-to-near infrared colors such as R — K > 
5 — 7, I — K > 4 — 6 (Vega magnitudes). Spectroscopic studies have shown that ob- 
served EROs can be classified as either dusty star-forming galaxies, old passive galaxies or 



AGN flSpinrad et al.lll997l : iDev et al.lll999l : ICimatti et al. 



1999; Soifer et al. 1999 



Liu et al. 
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20001 . e.g.). IPozzetti &: Mannuccil (120001 ) have shown that the old stellar population is ex- 
pected to have a strong spectral break between 3000 and 5000 A, while the effect of dust 
reddening in a younger population would produce a smoother spectra. This effectively 
segregates the ERO population into dusty star-forming and passive old galaxies (hereafter 
dERO and pERO, respectively) at 1 < z < 2 in the RJK or UK colour diagram. Man- 
nucci et al. (2002) have shown that the relative abundances of dEROs and pEROs is 
roughly equal using a s ample of 57 EROs {K[ < 20 and R — K' > 5.3) based on the 



classification method of IPozzetti fc Mannuccil (120001 ) (Their sample was classified as 37% 
elliptical, 37% star-forming, 9% stars with 17% remaining unclassified). Subsequent spec- 
troscopy, morphological and X-ray observations have shown that the comoving space densi- 
ties of elliptica 



(ICimatti et al 



-like and disk-like EROs ar e comparable over the magnitude range Ks < 20 
Moustakas et al.ll2004f) and that only a small fraction (5%-1 0%) of EROs 



2002 



host AGN (IRoche et all I2003T) . The additional properties of strong clustering (IDaddi et al. 



2000l ; McCarthy et al. 120011) and the relation of dEROs to (ultra) luminous IR-galaxies has 
also been examined (Elbaz fc Cesarskyl 12003). Indeed dEROs have been detected in the 



i nfrared with both t he ISO observatory (jElbaz et al.ll2002l ) and the Spitzer space telescope 
(jWilson et all l2004j ) .with more than 12% of the 8-24/zm sources being classified as dEROs 
in Spitzer observations of the lockman Hole. dEROs have also been detected at longer 
submillimetre wavelengths by the Submillimetre C ommon User Bolometer A rray (SCUBA) 
on the James Clerk Maxwell Telescope in Hawaii (IFrayer et al.ll2003l . 12004 ) . If dEROs are 
detected in significant numbers at infrared-submillimetre wavelengths then they may also be 
significant contributors to the dusty cosmic star-formation rate in the Universe. 



Recently, IDaddi et all (120041 ) have suggested a new criterion to isolate star-forming and 
passive evolving galaxies at 1.4 < z < 2.5 in the color-color parameter space based on B-, z- 
and K -band photometry. Previous studies of these BzK selected sources have shown that 
the star-forming galaxies (hereafter sBzKs) have star formation rates of ~ 1OOM yr _1 and 
stellar masses of 10 10_11 M Q , whilst the passive evolving galaxies (hereafter pBzKs) are found 



to ha ve stellar masses of > 1O U M (IDaddi et alJl2004 l2005t iKong et al.ll2006l: lHavashi et al 



2007 ). Both pBzK s and sBzKs are also found to be strongly clustered (IKong et al 



Hayashi et al 



2006 



2007J) 



Distant Re d Galaxies ( D RG) are identified by their observed red near-infrared colors, 
J — K > 2.3. iFranx et al. ( 2003 ) argue that the DRG color criterion is designed to be 
sensitive to those galaxies with a large 4000 A break and produces a sample that is mainly 
populated by galaxies at z > 2. Subsequent analysis has concluded that DRGs have high star 
formation rates of 100 — 400 M m yr" 1 , high stellar masses of (1 — 5) x 10 n M^ and are actively 
forming stars at z ~ 2 — 3.5 (IRubin et al.ll2004l ; Ivan Dokkum et al.ll2004l ; iKnudsen et al.ll2005 
Papovich et al.ll2006l . e.g.). 



-4- 



The contamination by galaxies at lower redshifts of the high-z populations selected 



via the DRG colour c riterion has also been discussed in recent papers (IGrazian et al.l 12006 



Papovich et al.l 120061 ) and it has been shown that the high-z DRG population are more 
luminous and exhibit stronger clustering than that of the low-z population that is also 
selected with the DRG colour criterion. From these studies, these high redshift galaxy 
populations are inferred to be the progenitors of the present-day early type galaxies. 

In this work, we focus on the number counts of such K- band selected high redshift 
galaxy populations. Based on a very deep sample (Kg < 22) of 198 EROs from the 
ESO/GOODS data in the CDF-S, iRoche et al.l (|2003h show that the ERO number counts 
exhibit a change of slope flattening at K > 19.5 from d(log N) /dm = 0.59 to 0.16. Re- 
cent results of a deep survey with the Subaru/MOIRCS instrument in the GOODS-N field, 



Kaiisawa et al 



(IFoucaud et al. 



200 



6; 



(20071 ) shows that the nu mber counts of DRGs also turn over at K ~ 22 



Grazian et al.ll2006l see also). Although the number counts of sBzKs 



increase in number sh arply toward fainter magnitudes, th ose of the pBzKs also clearly turn 
over around K ~ 19 (IKong et al.l 120061 ; lLane et al.l 120071 ). This indicates that these num- 
ber counts are directly sampling the luminosity function (LF) at a specific epoch. In this 
scenario, the corresponding absolute magnitude at the flattening of the number counts will 
be close to the value of determined for the K -band LF oy er the same redshift range 
( Kashikawa et al. 2003 ; Pozzetti et al. 2003 ; Saracco et aL 12006). Furthermore the trend of 
a late-type galaxies to poses s a steeper slope of t he LF (IFolkes et al.lll999l : iMadgwick et al. 



fr 

20021 ; iNakamura et al.l 120031 ; iDahlen et al.l 120051 ) is also consistent with that of the sBzK 



counts exhibiting steeper slope at the faint end. 

Note that some studies have claimed that the ERO counts are lo wer than those pre- 
dicted by pure luminosity evolution (PLE) models. IRoche et al.l (120021 ) derived the number 
counts of EROs using a sample (K < 21) of 158 EROs in the ELAIS N2 field and re- 
ported lower counts than models in which all E/S0 galaxies evolve according to PLE. On 
the other hand, conflicting r esults have also been found from brighter samples of EROs. 
Vaisanen &; Johansson! (120041 ) show that the cum ulative number counts of ERO s (K < 17.5) 
over a 2850 arcmin 2 area within the ELAIS field (IRowan-Robinson et al.ll2004l ) are well fited 
by PLE models. However, several results implying a lower density of early-type galaxies at 
z > 1 have been reported from morphological studies. For instance, iKajisawa &; Yamada 
(120011 ) reported that the comoving density of My < —20 galaxies decreases at z > 1 in each 
morphological category based on the results of the HST WFPC2/NICMOS archival data 
of the HDF-N. In particular, they showed that the number density of the bulge-dominated 
galaxies conspicuously decreases to 1/5 of the local value between the 0.5 < z < 1.0 and 
1.0 < z < 1.5 bins. Other morphological studies have also revealed an apparent lack of mas- 
sive ellipticals at redshifts greater than z ~ 1.5 (IFranceschini et al.lll998t iRodighiero et al. 
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20011 . e.g. 



Hierarchical models of galaxy formation predict a significant decline in the comov- 
i ng density of ellipticals with redshift, as they should form through merging at z < 2 



(W hite & Rees 



1999 



ICole et al. 



978l:IWhite fc Frenklll99ll ; lKauffmann et al.lll996l : lBaugh et al. 1996; S omerville fc Primacl 



2000l . e.g.) and a measure of a decline in their comoving density would pro- 



vide valuable proof of these models. It is therefore important to study the counts of these K 
-band selected high redshift galaxy populations because these samples will provide us with 
the clues for the formation and evolution of galaxies. 

We describe our observations in the AKARI north ecliptic pole (NEP) field in §[2j 
We then present our sample of the K -band selected high redshift galaxy populations in 
§|HJ These observed number counts are presented in §HJ In §0 we discuss the relationship 
and trend of the ERO, BzK and DRG count slopes, the results of comparing with galaxy 
count models and implications from other observational studies. Finally, the summary of 
this research is given in §|6j Throughout this paper a concordance cosmology of Ho = 70 
kms _1 Mpc _1 , fl m = 0.3 and Q\ = 0.7 is assumed in a greement with the results for Typ e 



la supernovae (IRiess et al.lll998l : IPerlmutter et al.lll999l ) and WMAP flSpergel et all 120031 ). 
Unless otherwise stated, we use the Vega magnitude system. 



2. AKARI NEP field 



The Japan A erospace Explora tion Agency (JAXA) launched the infrared astronomy 
satellite AKARI flMurakamil booi formerly ASTRO-F) on February 22nd 2006 (Japan 
Standard Time). 

Due to the fact that the orbit of AKARI is Sun-synchronous polar, AKARI can often 
obtain deep exposures at the e cliptic poles. Such a de ep survey program is currently under 
way in the AKARI NEP field (IMatsuhara et al.l 120061 ). A portion of the field has also been 
covered by deep B, V, R, i' and z' optical observations with Subaru/Suprime-Cam, reaching 
B = 28.4, R = 27 A and z' = 26.2 AB magnitude in 5 a (Wada et al. 2007 inpreparation) 
and J, Ks near-infrared observations with KPN02.1 m/FMAMINGOS, where t he deepest 
regio ns reach to J = 21.85 and Ks = 20.15 Vega magnitude in 3 a, respectively (jlmai et al. 



20071 ). In this paper, since we focus on the K s -band selected high redshift galaxy popula- 
tions, we adopt the catalog of the region reaching to K s ~ 20 Vega magnitude in 3 a. The 
total effective area is 540 arcmin 2 . 



- 6- 



3. Sample Selection of EROs, BzKs and DRGs 



3.1. The ERO Sample 

The selection of EROs depends on the depth of the photometry, the available filters and 
assumed color criteria. In this paper, alth ough EROs are define d as ob jects with R — K s > 
5.3 following the classification scheme of iPozzetti fc Mannuccil ( 120001 ) . we also present an 
alternative ERO selection criteria of R—Ks > 5.0 to compare our results with other published 
ERO counts. 

Figure [1] shows the R — Ks color versus Ks magnitude for all detected objects in the 
AKARI NEP field. In total, we identify 308 (442) EROs with R-K s > 5.3 (R-K s > 5) to 
Ks = 19.0. We also calculated completeness correction factors as a function of Ks magnitude 
by applying Monte Carlo simulations to the real images (see llmai et al.l (120071 ) for details 
of the method). Above the magnitude of K s = 19 the catalogue appears more than 70% 
complete. Though the ERO candidates include the R -band undetected sources, this is less 
than 4% of the total ERO counts. Some faint galaxies may be identified as point-like objects 
because their extended structure vanishes into the background noise. The cross symbols in 
Figure [T] represent point-like objects which are defined as less than 5.4 pixels (corresponding 
to 1.08") in the z' -b and image (a verification of the star and galaxy separation is described 



in 



Imai et al.l (120071 ) ) . Our result shows that two point-like objects in the magnitude range 
17 < Kg < 19.0 satisfy the ERO selection criteria of R — Ks > 5.3. Although, we do not 
include these objects in the ERO counts, their contribution is ~ 0.6% of the total ERO 
counts and therefore will not sig n ifican tly affect the results. Figure [2] shows the separation 
method of IPozzetti &: Mannuccil (120001 ) . The extracted numbers of pERO and dERO are 
167 and 141 to K s = 19.0, respectively. 

Note that an alternative selection method to segregate the pERO and dERO popu- 
lations would be a direct measurement of their mid- & far-infrared flux. Using Spitzer 



Stern et al.l (120061 ) observed the two archetypal extremely red objects HR10 (dERO) and 
LBDS 53W091 (pERO) from the IRAC 3.6 ^m band through to the MIPS 160 ^m band. 
Although both sources were well detected in the four shorter IRAC bands, only HR10 (i.e. 
a dusty star-forming galaxy) was detected in the MIPS 24 & 70 /im bands. Obviously, mid- 
far-infrared detections of EROs provides a powerful alternative segregation criteria for dusty 
EROs and we will present similar AKARI mid-infrared fluxes for our sample in a later work 
on completion of the AKARI NEP survey. 
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3.2. The BzK Sample 



We select sBzKs and pBzKs by using the color criteria of iDaddi et al.l (12004 ) . Figure [3] 
shows the BzK color diagram of Ks -band selected objects in the AKARI NEP field . In 
order to apply the BzK selection criteria in a way consistent with IDaddi et al.l (120041). w e 
normalize the stellar sequence in the AKARI NEP survey to that of IDaddi et al.l (120041 ). 
resulting in a shift of 0.1 mag to fainter Ks magnitudes. From the results, we further selected 
38 p BzKs and 137 sBz Ks to Ks = 19.0, above which the catalog is more than 70% complete 



sec 



Imai et al.l (120071 ) for details). Seven point-like objects also satisfy the sBzK criteria. 
Since their contribution is less than 5% of the total sBzK counts, we do not include them in 
the sBzK counts. 



3.3. The DRG Sample 



We selected the D RG sample using the J — Ks > 2.3 criterion (IFranx et al.l 120031 ; 
van Dokkum et al.ll2003l ). Figure H] shows the J — Ks color against Kg magnitude for all 



detected objects in the AKARI NEP field. The objects shown by the cross symbols are 
point-like objects. For the DRG candidates, there is no contamination from point-like objects 
using the J — Kg > 2.3 criterion. However the J -band images are not deep enough to find 
counterparts for all the Ks -band selected objects. We therefore restrict the DRG sample to 
be brighter than Ks = 18.6. T his threshold inc l udes a ll DRGs bluer than J —Ks ~ 2.9 which 
is almost the same criterion of iFoucaud et al.l (120061) and provi des 64 DRGs to K s = 18.6, 
above which the catalog is 90% complete (see llmai et al.l (120071 ) for details). 



4. EROs, BzKs and DRGs counts 



In this section, we present the Ks -band number counts for the ERO, BzK and DRG 
populations (the relationship between th counts are discussed in §5.ip . The raw numbers of 
these K s -band selected high redshift galaxy populations and completeness corrected counts 
per square degree per magnitude are listed in Table HI [2, [3] and HI 

The ERO counts are presented in Figure [U which shows the results from the different 
selection criteria along with a compilation of published ERO counts. The E RO counts in the 
AKARI NEP field are in good agreement with those of IDaddi et al.l (120001 ) where the EROs 
were selected using an identical criteria. Howev er our numbe r coun ts in the faint magn itude 
range Ks > 18.5 are slightly lower than those of iRoche et al.l (120031 ) ; iKong et al.l (120061 ) who 
utilized different filter bands to select the EROs. 
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Although there are many examples of t otal ERO counts and the segregation of the 



two constituent populations in the literature ( Roche et al 



Vaisanen &: Johansson! I2004J : iKong et al.l 12006 



2002 



20031 : iMivazaki et al.l 12003 



Simpson et al.ll2006l ). the individually segre- 



gated source counts of both pERO & dERO populations have yet to be coherently presented. 
In this work, we also present the observed counts of pEROs and dEROs separately in Figure 
|6j In the bright magnitude range, the pERO counts are high er than that of the dEROs; by 
a factor of ~ 1.5 at Kg ~ 17.6. This result is consistent with IVaisanen fc Johansson! (120041 ) 
who concluded that the fraction of pEROs range from 65-80% at K < 17. 

Figure [7] and [8] shows th e number counts of p BzKs and sBzKs in the AKARI NEP field 



together with the results of iReddy et al.l (120051 ) ; IKong et al 



Our pBzK counts are 



broadly consistent with the co unts in the Daddi field of IKong et al.l (120061 ) while the number 
counts in the Deep3a field of IKong et all (bood ) and UDS EDR field of Lane et al.l bpOTh 
are higher by a factor of 3 at Ks ~ 18.5. 

The discrepancies between the field to field counts could in some part be attributed 
to varying selection criteria, however another possibility is variation in the number den- 
sities due to the effects of cosmic variance. For small size fields and strongly clustered 
populati ons such as EROs, th e uncertainty due to the cosmic variance can be as high as 
30-40% (ISomerville et al.l 120041 ). where as larger survey areas such as the UDS EDR field of 



Lane et al 



3-^ 



(120071 ) of 0.6 deg 2 are le ss susceptible with an estimated cosmic variance at the 
o level over an interval 5z ~0.2 (IPeacock. Doddslll994l ). 



Combining the data with lReddv et al.l (120051 ) . IKong et al.l (120061 ) and lLane et al.l (120071 ). 
the slope of the pBzK counts shows a turn over at Kq ~ 19.0. On the other hand, our sBzK 
counts are in good agreement with IKong et al.l (120061 ) in both the Deep3a and Daddi fields 
at all magnitudes 17.5 < Ks < 18.8 within the error bars. The sBzK counts are more than 
5.8 times higher than the pBzK counts at Ks ~ 18.6 and show a steeper slope than those of 
the pBzKs in the faint magnitude range. 



Fi nally we present the D RG counts in Figure [9] together with the results of iReddv et al. 



( 2005 ); iFoucaud et al.l (120061 ). The DRG counts are in good agreement with iFoucaud et al. 
(120061 ) and the slope of the counts do not show a turnover to K s < 18.5. 



5. Discussion 

We have presented ERO, BzK, DRG counts in §H In this section, we discuss the 
relationship between, and the trend of the ERO, BzK and DRG count slopes (I5.ip . the 
results of comparing our counts with evolutionary models for the passive early-type galaxy 



- 9- 



population (I5.2p and the dERO star-forming population (15.31) . 



5.1. Difference between the ERO, BzK and DRG counts 



As previously indicated by iRoche et al.l (120031 ) in other ERO fields, the ERO counts in 
the AKARI NEP field (Figure^]) also clearly show a slope change at K$ ~ 18.0 below which 
the counts exhibit a shallower slope. On the other hand, the DRG counts show no sign of 
any flattening of the slope at K$ < 18.5, although it should be noted that the DRG counts 
in the AKARI NEP field are drawn from limited data across only three magnitude bins. 

Ho wever a shallower slope has been reported at fainter magnitudes f rom t he deeper DRG 



counts (IFoucaud et al.l 120061 ; iGrazian et al.l 120061 ; IfCajisawa et al.l 120071 ) and iFoucaud et al. 
(120061 ) suggest a break feature in the slope at Kab ~ 20.5. Although EROs and DRGs both 



include dusty star-forming and passive old populations, these results suggest that the red 
galaxies contribute less significantly to the galaxy counts in the fainter magnitude range. 

The pERO counts show a turn over at K$ ~ 18.5 whilst the pBzK counts show a turn 
over at Ks ~ 19.0, shifting to slightly fainter magnitudes. This is reasonable since the pBzK 
criterion is generally sampling a higher redshift domain than the pERO criterion. Since the 
redshift distributions of the pERO and pBzK counts are well-defined, these populations are 
independently selected passive old galaxies in different redshift slices. Therefore the pERO 
and pBzK criteria can provide a good opportunity to investigate the evolution of early-type 
galaxies. We will investigate the difference between these counts of passive old samples and 
produce a number count model for the early-type galaxies in next subsection. 



5.2. Comparison with early-type galaxy count models 



In this section, we construct a number count model for the early-type galaxies in order 
to explain the pERO and pBzK coun ts. We primarily adopt the number count model based 
on the formulation of I Gardner! (119981 ). This model incorporates the evolution of the LF for 
each galaxy population within a backward evolution framework. Therefore, we may control 
the evolution of each galaxy type using the type- dependent LF. 



The number count model depends mainly on the three parameters of the ISchechter 

<j>(L)dL = f(±-) exp 



fll976j ) LF 



L\dL 
L*) 



1 
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with 



L 



10' 



0.4(Af*-Af) 



(2) 



The parameter M* is the characteristic magnitude and the parameter a is the faint-end slope 
of the LF. Although the parameter 0* is a normalization coefficient with dimensions of the 
number density of galaxies, we re-define <fr* as follows, in order to include a parameterization 
of evolution in the number density, 



0*(O)(l + z) 



(3) 



where m is a free parameter controlling the number density evolution. In the case of m < 0, 
the number density is decreasing to higher redshift (hereafter 'negative density evolution'). 



In the past decade, there have been a number of determinations of t 



galaxies based on both ph otometric (IBalogh et al.ll200ll: iKashikawa et al 



2005 



2003 



Saracco et al. 



Pozzetti et al. 



2006) and spectroscopic redshifts (IGlazebrook et al. 



he NIR local LF of 
20031 ; ICaputi et al. 



1995 



Huang et al. 



20031 ). From optical surveys, it has also been shown that each galaxy 



type has a characteristic LF shape (IFolkes et al.lll999l ; lMadgwick et al.ll2002l ; iNakamura et al. 



20031 ; iDahlen et al.ll2005l ). However LFs determined with optical surveys may be affected by 
star-formation activity and dust extinction, 



we refer to the type-dependent i^-band luminosity function (LF) of iKochanek et al. 



( 20011) based on th e 2MASS Second Inc r ement al Release Catalog of Extended Sources 
( Jarrett et al. _ 200(3) . H owever Cole et al. ( 2001 ) point out that the i^-band photometry 



of lLoveday et al.l (120001 ) has better signal-to-noise and resolution than the 2MASS images 
and enables more accu r ate 2M ASS magnitudes to be measured than the original default 
magnitudes iFrith et al.l (120051 . see also). We therefore normalized the characteristic mag- 
nitude Mt of the number count model at K < 13.25 of the 2M ASS counts measu red by 
Frith et al.l (120051 ) using a similar magnitude estimator to that of ICole et al.l (j200ll ). The 
three parameters for the LF of early-type galaxies are shown in Table [5j 



We use the spectrophotometric population synthesis models of lBruzual fc Charlotl (120031 ) 
with a Salpeter stellar initial mass function to compute the evolution of the galaxy luminosity 
as a function of the age, which we refer to as our baseline Pure Luminosity Evolution (PLE) 
case. This theoretical spectral energy distribution (SED) model is ideal for the approach to 
number count modelling within the backward evolution framework because it returns the 
SED of a galaxy with any star formation history as a function of cosmic time for a given 
formation epoch. The basic characteristics of the SED of early-type galaxies i.e., metalicity 
and star formation rate, assuming pure luminosity evolution, are also listed in Table 



The galaxy formation epoch is an important model parameter, as any delay in the galaxy 
formation epoch will shift the count slope towards brighter magnitudes and lower number 
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densit ies. It is, however, poorly understood and rather loosely constrained. iHolden et al. 
( 120041 ) reported a lower limit for the elliptical galaxy formation epoch from a study of the 
color-magnitude relation of elliptical and SO galaxies in six clusters obtained using imaging 
on the HST. They reported the oldest stars in the elliptical galaxies appear to have formed at 
z > 3. We therefore calculate our number count models for the early-type galaxy population 
for three different formation epochs at Zf orm = 4, 5 and 6. 

We first compare the calculated early-type number counts with the observed pERO 
counts as shown in Fgure dUJ The effect of the strong negative number density evolution 
is seen as a shift of the counts to lower number densities with magnitude. In the bright 
magnitude range Ks < 17.0, since the pEROs sample is small in number, the simple PLE 
model without density evolution (corresponding to m = 0) cannot be ruled out. However 
in the fainter magnitude range Ks > 17.0, the simple PLE model is clearly inconsistent 
with the pERO counts: the faintest pERO counts at Ks ~ 19.1 are more than 6.7 times 
lower than the PLE models with Zf orm = 5. Note the pERO counts are higher than the 
negative density evolution model with m = —1.4, this is because the brighter samples of 
pEROs contain many low redshift galaxies (jVaisanen fc Johansson! 120041 . e.g.). We conclude 
that the pERO counts favor the negative density evolution models. 

We also compare the early-type number count model with the observed pBzK counts. 
The simple PLE models with m = exceed the observed counts at all magnitudes while the 
negative density evolution models with m = —1.4 agree w ell with the bright p BzK counts 



i n the AKARI NEP field and the deeper pBzK counts of iReddy et al.l (120051 ); iKong et al. 

(hood ). 



The number count model with negative number density evolution (m = —1.4) shows 
better fits for both the pERO and pBzK counts. Since independently selected passive old 
samples both show signs of negative density evolution, these results indicate that the overall 
early-type galaxy population may exhibit negative density evolution. 

Indeed negative density evolution has also been reported from studies of the evolutio n 



of the K -band LF at various redshifts: z < 1 (Glazebrook et al 



z < 2 (Bolzone 



(ISaracco et al. 



1995 



Feu 



la et alil200~2l ; |Pozzetti et al.ll2003h . z < 3.5 jKashikawa et al 



ner et all 120031) 



2003|) and z < 4 



2006h . We compile these results in Figure [12], and compare wit h the negative 



densit y evolution model. The best fit is obtained for m ~ —1.3. Note that iPozzetti et al. 



(120031 ) report a contrary behavior, however, their large errors due to the limited sample size 
do not rule out a decrease in the comoving number density with redshift. To summarize, our 
counts are consistent with strong negative density evolution (m = —1.4) in the early-type 
galaxy population color selected (pERO and pBzK) sample, whiles the total rest-frame K 
-band LF also exhibits negative density evolution with redshift with a best fit parameter 
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m ~ 



-1.3. 



5.3. Comparison with dusty star-forming galaxy count models 



In this section, we will construct a simple number count model for dusty star-forming 
galaxiesto represent the dERO counts. We assume that the nature of the dERO's are high 
redshift starbursts with their redness being due to the extreme dust extinction within the 
host galaxies. It follows that if indeed the dEROs are significant sites of star-formation then 
these sources shou ld also exhibit significant emission at infrared-submillimetre wavelengths 
felbaz et al.ll2002h . 



In recent years a great deal of effort has gone into estimations of the fraction of dEROs 
making up the far-infrared-submillimetre counts, however these estimates have been ham- 
pered by differing selection criteria. Estimates of the fraction of dEROs in the bright 
(S 850 ^ m > 2mJy) submillimeter population hav e been made in th e SCUBA -CUDSS fields 
where 10-20% of sources were identified as EROs ( Webb et al. 12003 ). although Clements et a! 
(120041 ) found a higher fraction of 8/14 in the CUDSS 3hr field sources could be classi- 
fied as EROs. In the SCU BA-8mJy fields, 6 /30 sources (6/1 7 with optical identi fications) 
were identified as dEROs (jlvison et al.ll2002l ). More recently, iTakagi et al.l (120071 ) have es- 
timated a dERO con tribution of 10±7% (for Kg < 20.3) in the SCUBA-SHADES survey 
(IMortier et al. 1120051 ) although this fraction may increase to fainter K magnitudes. However 
within the bright submillimetre fields, the average source density of near-infrared selected 
(p+d) EROs is ~18,000 per square degree, i.e. significantly higher than the source den- 
sity of bright submillimtre galaxies (^thousands per square degree). Therefore, for the 
dEROs to be explained by the submillimeter source population there must exist a larger 
contribution at fainter submillimetre fluxes, and indeed, deeper submillimetre observations 
{.S^OfMm <2 mJy) have implied a much higher fraction of EROs. Using the HDF-N SCUBA 



super-map, iPope et al.l (120051 ) discovered a much higher contribution of EROs than in the 
brighter submillimetre surveys and moreover classified all submillimetre detected EROs as 
dEROs. At submillimetre fluxes of Sg^m ~lmJy, the source density of submillimetre 
sources becomes sufficiently high to account for a ll EROs. To investiga te the contribution 
of EROs to the fainter submillimetre population, iKnudsen et al.l (120051 ) performed submil- 
limetre stacking analysis of red galaxies around distant clusters finding an average flux of 
Ssso/im ~lmJy and an estimated contribution to the total flux from faint sources between 
0.5 < Ss5o^m < 5m Jy of 50%. Similar studies have yielded consistent results with average 
fluxes of 1.58±0.1 3mJy and contributi ons to the emission over the 0.5 < Sssotim < 5mJy 
flux range of 50% (IWehner et al. 1120021 ). and 0.4±0.07mJy, with a contribution to the total 
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emission over the 0.5 < Ss^o^m < 5mJy range of 10% (IWebb et al. 112004 ) respectively. The 
differences are most likely due to the slightly different ERO depth and selection criteria and 
field to field variations due to the strong clustering nature of dEROs. 

Under the above assumptions, we model the K-band dusty star-forming ERO popu- 
lation as a fraction of the dusty far- infrared - submillimetre po pulation. We assume the 



submillimetre luminosity function of ISerjeant fc Harrison! (120051) that w as derived from a 
combination of the I RAS point source c atalogue (ISaunders et al.l 120001 ) and the SCUBA 
Local Galaxy Survey (IDunne et al.l 120001 ). The luminosity function is a double power law 
(somewhat broader than the traditional Schecter function in the higher luminosity regime) 
and is defined as; 



(f>(L)dL 



dL 



(4) 



where L* is the characteristic luminosity (with corresponding characteristic magnitude, 
M*, given by L/ L*=IQ 0A ( M *~ M ^ and (p* the number density normalization coefficient at 
4>(L*). The parameters, (3 & 7 control the faint and bright end slopes of the luminosity 
function. The adopted luminosity function parameters are shown in Table [6j We define 
three bolometric luminosity regimes for our luminosity function of L/L & >10 10 , 10 11 , 10 12 
corresponding to star- forming infrared galaxies (STFG), luminous infrared galaxies (LIRG) 
& ultra- luminous infrared galaxies ( ULIRGs) r espec tively. These populations are modelled 
within the evolutionary framework of iPearsonl (120081 ) which includes both luminosity evolu- 
tion (k) and density evolution (g) as double power laws, of the form (1 + z) k,g , rising steeply 
to a redshift of approximately one and shallower thereafter. The evolutionary strength and 
form is defined by a pair variables for both the luminosity and the density evolution: the 
power index from redshift zero to redshift unity, k\,g\ and an additional power index, fc2,5 l 2, 
to a higher reds hift Z2 ~ 3 .0, af ter which the evolutionary power is assumed to decrease ex- 
ponentially (see IPearsonl (120081 ) for explicit details of the model). The parameterization of 
the power laws are given in Table [7J 

Our source SEDs are selected from the spectral libraries of Takagi et al. ( 2003al lbl) which 
provide excell ent fits to the SEDs of submillimetre galaxies from submillimetre to near-infared 
wavelengths (ITakagi et al.l I2004J ) . These models assume a Salpeter IMF of index 1.35 from 
O.1-6OM and parameters for the star formation timescale (O.lGyr), metallicity (assumed to 
be 0.1), the compactness of the starburst region (O) and the starburst age. The variation in 
the starburst SED is explained by the difference in the starburst age and the compactness 
of the starburst region. The starburst age ranges from 0.01-0.6Gyr. The compactness factor 
O ranges from 0.3 - 5.0 and can be considered as a measure of the optical depth. In order to 
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select SEDs appropriate for our modelling of the dERO population we identify those SEDs 
from within the spectral library that specifically have red colours R — K >5.3 in the redshift 
range 1-2. The selected SEDs for our three populations are tabulated in Table Note that 
for our L/Lq >10 12 population (i.e. ULIRGs) we adopt t he model fit to the archetypical 
dERO HR10 (ERO J164502+4626.4. Illu fc Ridgwavl Jl994h l which has been shown to be a 
distant Arp220 type ULIRG analogue felbaz et al.ll2002h . 

Using the above luminosity function and SEDs we model the observed counts of dEROs 
in the AKARI NEP field by assuming ERO fractions of the submillimetre population of 15, 
25 & 35 percent respectively. The resulting counts are shown in Figure [131 The i^g-band 
dERO number counts are well fitted by a submillimetre evolutionary model assuming that 
around 25 percent of the sub-millimetre galaxy population to faint fluxes are dEROs. The 35 
percent model over-predicts the dERO counts at K$ magnitudes fainter than 17, whilst the 
15 percent model cannot account for the fainter dERO source counts at K$ >18.5. In Figure 
[Hwe also plot the individual components (STFG L/L >10 10 , LIRG L/L >10 n , & ULIRG 
L/L Q >10 12 , see Table [6]) corresponding to our best fit 25 percent model given in Figure [TBI 
In our best fit 25 percent model, the bulk of the brighter (K$ <19) counts originate from 
the L/Lq <10 12 population, i.e. that HR10 type dEROs do not contribute significantly to 
the dERO source counts . Not e that this result is consisten t with the sub-mm stacking of 
ULIRGs of iMohan et al.l (120021 ) and the work of ISmail et al.l (119991 ) who concluded that the 
contribution of HRlO-like dEROs to the submillimetre population was constrained to 10-20 
percent. The brighter counts at (K$ <18) are composed of LIRG sources with luminosities 
10 11 < L/L Q <10 12 along with a non-negliable fraction of weaker L/<10 11 L Q star-forming 
galaxies. We also predict that by K$ >19.5, the dERO source counts are dominated by the 
faint submillimetre population of moderate star-forming galaxies. Note that if this is the 
case, then the fainter dEROs are examples of less luminous systems responsible for a large 
fraction of the star- formation rate at high redshift. We also predict that, corresponding to 
the decline of the brighter systems, the counts of dEROs should turn over in the magnitude 
range 7^=19-20. 



6. Summary 

In order to seek clues to understand the origin of the change in the K -band count slope 
and the evolution of high redshift galaxy populations, we have investigated the number counts 
of K -band selected high redshift galaxy populations: EROs, BzKs and DRGs. Moreover 
we have further segregated these populations into their passive elliptical and dusty evolving 
star-forming constiuents. Our catalogue in the AKARI NEP field contains 167 pEROs and 
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141 dEROs (K s < 19.0), 38 pBzKs and 137 sBzKs (K s < 19.0), 64 DRGs (K s < 18.6). From 
these large samples, we confirmed that the number counts of both the pEROs and pBzK 
have shallower slopes at fainter magnitudes and turnovers at Kg ~ 18.5 and Ks ~ 19.0, 
respectively. On the other hand, the sBz K counts show a ste e per slope whic h dom inates in 



the fainter magnitude range as shown by iKong et al.l (120061 ); lHayashi et al.l (120071 ) . These 



results strongly support that the passively evolving population contributes less to the galaxy 
counts at fainter magnitudes. 

For the passively evolving pERO and pBzKs populations, from a comparison with a 
backward evolution model allowing number density evolution, we conclude that the PLE 
model without density evolution is at odds to explain the observed counts at fainter magni- 
tudes, whilst the model including negative density evolution is consistent with the observed 
counts of the passively evolving populations. 

We have also modelled the dERO population with a simple backward evolutionary 
model assuming that the dERO population can be explained by the source densities of 
the far-infrared - submillimetre populations. Our best ft to the dERO source counts is 
found by assuming a 25 percent contribution from submillimetre star-forming galaxies with 
the majority of brighter Ks -band detected dEROs having luminous, as opposed to ultra- 
luminous counterparts. With such a significant contribution by dEROs to the fraction of 
dusty galaxies at far-infrared-submillimetre wavelengths, we may expect a correspondingly 
significant contribution, from this population, to the cosmic star-fromation history of the 
Universe. 

Further investigation of these high redshift colour selected populations (ERO, BzK, 
DRG) will not only produce constraints on the formation epoch and evolution of giant 
galaxy systems but also provide a bridge to connect the optical and infrared star-formation 
history of the Universe. A discussion of the AKARI mid-infrared emission from our sample 
will be given in a forthcoming work. 
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Fig. 1. — R — Ks color against K$ magnitude in the AKARI NEP field. The open boxes and 
filled circles represent dusty star-forming and passive old EROs selected with the R — Ks > 
5.3 criterion (the horizontal dashed line). The classification of dusty star-forming and passive 
old EROs is given in Figure [2j The dots are other extended objects while the cross symbols 
are point-like objects. The vertical solid line represents the magnitude limit at Ks = 19.0. 
The thin-dashed line shows the upper bound of the R—Ks color due to the survey magnitude 
limit; galaxies above this line are too faint to be detected in the R -band (R > 27.0 mag). 
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Fig. 2.— J — K s vs. R - K s colors of EROs in the AKARI NEP field. The dusty star- 
forming ERO / passive old ERO separation is made based on the color discriminator from 
Pozzetti fc Mannuccil (j2000l ). The open squares and filled circles represent dusty star- forming 
and passive old EROs, respectively. 
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Fig. 3. — B — z' vs. z' — Ks colors of Kg -band s elected objects in th e AKARI NEP field. 
The discriminating color criterion is adopted from lDaddi et al.l (j2004l ). The open boxes and 
filled circles represent the star-forming and passive old BzKs, respectively. The dots are other 
extended galaxies while the cross symbols are point-like objects. The dots and two dotted 
lines represent the (V — Ks)ab = 2.5 and [z! — Ks)ab — (B — z')ab = —0.2, respectively. 
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Fig. 4. — J — Kg colors against Kg magnitude in the AKARI NEP field. Filled circles 
represent DRGs selected with the J — Ks > 2.3 criterion while dots are other extended 
galaxies. The cross symbols are point-like objects. The dashed and solid lines represent the 
J — Ks > 2.3 criterion and the magnitude limits at Ks = 18.6, respectively. The dotted line 
shows the upper limit of J — Ks color due to the J -band limit; galaxies above this line are 
too faint to be detected at J -band (J > 21.5 mag). 
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Fig. 5. — Differential galaxy number counts for EROs. The filled circles and diamonds are 
the counts in the AKARI NEP field for EROs extracted according to the different selection 
criteria R — K$ > 5.0 and R — K$ > 5.3, respectively. Error bars are calculated from 
the number of galaxies in each magnitude bin using Poisson ^/N counting stati stics. Other 



symb ols in the figure correspond to the coun ts in the Deep3a a nd Daddi field (IKong et al. 



20061 ). CDF-S (IRoche et al.l2003h . ELAIS N2 (IRoche et al.ll2002h . The open circles represent 
the results of baddi et all (l200oh . 
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Fig. 6. — Differential galaxy number counts for dusty star-forming and passive old EROs 



i n the AKARI NEP field based on the separation criteria given by iPozzetti fc Mannucci 



( 120001 ) . The open squares and filled circles correspond to dusty star-forming and passive old 
ERO counts, respectively. Error bars are calculated from the number of galaxies in each 
magnitude bin using Poisson \fN counting statistics. 
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Fig. 7. — Differential galaxy number counts for passive old BzKs. The filled circles corre- 
spond to the counts in the AKARI NEP field. Error bars are calculated from the number 
of galaxies in e ach magnitude b in using Poisson yN c ounting statistics . The counts in 
the UDS EDR Lane et al.l l2007h. Deep3a, Daddi field fcong et al.l bo06h and GOODS-N 
jReddv et alJboOfih :ields are also shown with other symbols. 
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Fig. 8. — Differential galaxy number counts for star-forming BzKs. The filled circles corre- 
spond to the counts in the AKARI NEP field. Error bars are calculated from the number 
of galaxies in e ach magnitude b in using Poisson yN c ounting statistics . The counts in 
the UDS EDR Lane et al.l l2007h. Deep3a, Daddi field fcong et al.l bo06h and GOODS-N 
jReddv et al.lhoosh fields are also shown with other symbols. 
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Fig. 9. — Differential galaxy number counts for DRGs. The filled circles are DRG counts 
in the AKARI NEP field. Error bars are calculated from the number of galaxies in each 
magnitude bin using Poisson yN countin g statistics. Other sy mbols in the figur e denotes the 
counts in UKIDSS-UDS, UKIDSS-DXS fIFoucaud et alJbooa i and GOODS-N faeddv et al. 
200fih fields. 
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Fig. 10. — K s -band number count predictions for the passive old EROs. Passive old EROs in 
the AKARI NEP field are shown as filled circles. The bold solid and dashed lines represent 
the PLE and the PLE plus density evolution model for a formation epoch Zf orm = 5.0. The 
accompanying upper and lower thin lines correspond to the different formation epochs of 
z — 3, 4, 5, respectively. 



31 




K s [mag] 

Fig. 11. — Ks -band number count predictions for the passive old BzKs. Symbols for 
observed counts are the same as those in Figure [71 The bold solid and dashed lines represent 
the PLE and the PLE plus density evolution model for a formation epoch Zf orm = 5.0. The 
accompanying upper and lower thin lines correspond to the different formation epochs of 
z — 3,4, 5, respectively. 
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Fig. 12. — Evolution of Kg- band luminosity function. The dotted and dashed lines represent 
the result of the negative density evolution given by eq. [3j m ~ —1.3 and -1.4, respectively 
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Fig. 13. — K s -band number count predictions for the dusty star-forming EROs. The lines 
represent different assumed sub-millimetre population fractions. The bold solid corresponds 
to a contribution of 25 percent of dEROs to the submillimetre population and the upper dot- 
ted and lower dashed lines represent contributions of 35 percent and 15 percent respectively. 
Dusty star-forming EROs in the AKARI NEP field are shown as open squares. 
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Fig. 14. — Kg -band number count predictions for the dusty star-forming EROs. The lines 
represent different assumed galaxy populations corresponding to the bolometric luminosity 
regimes of star-forming galaxies (STFG) L/L Q >10 10 (dashed), luminous infrared galaxies 
(LIRG) L/Lq >10 11 (dot-dashed), & ultra-luminous infrared galaxies (ULIRGs) L/L Q >10 12 
(dotted). The bold solid corresponds to the total contribution of dEROs (the same as in 
Figure [TBI) . Dusty star-forming EROs in the AKARI NEP field are shown as open squares. 
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Table 1: The number of ERO samples in the AKARI NEP field 



Ks bin center 


R 


-K s 


> 5 


R- 


- K s > 5.3 




N 


N 

1 y corr 


Error 


N 


N corr Error 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) (7) 


16.625 


47 


47 


18 


40 


40 16 


17.125 


110 


110 


30 


53 


53 19 


17.625 


350 


36 


50 


240 


250 40 


18.125 


790 


890 


80 


490 


550 60 


18.625 


1000 


1420 


120 


700 


990 100 



Note. — Effective area of NE and SE fields is 540 arcmin 2 . Unit of the counts is per square degree per 
magnitude. 

(2) and (5) Raw ERO counts. 

(3) and (6) Completeness corrected counts. 

(4) and (7) Poisson y/N statistics. 



Table 2: The number of sample of dusty star-forming and passive old EROs 
Ks bin center Passive old EROs Dusty star-forming EROs 

■Nraio -N corr Error -N ra ->M N corr Error 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


16.625 


27 


27 


13 


13 


13 


9 


17.125 


40 


40 


16 


13 


13 


9 


17.625 


150 


150 


30 


90 


100 


30 


18.125 


270 


300 


50 


210 


240 


40 


18.625 


310 


440 


60 


370 


520 


70 



Note. — Description of the columns are the same as in the footnote of Table [T] 
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Table 3: The number of sample of star-forming and passive old BzKs 



Ks bin center 


Passive old BzKs 


Star- 


•forming 


BzKs 




N 


N 

1 " corr 


Error 


N 


N 

1 " corr 


Error 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


17.125 








13 


13 


9 


17.625 


20 


21 


12 


70 


70 


20 


18.125 


50 


60 


20 


180 


200 


40 


18.625 


90 


120 


30 


490 


690 


80 



Note. — Description of the columns are the same as in the footnote of Table [P Numbers i n par enthesis 
include the point-like objects that also satisfy the BzK color selection criteria of , 



Daddi et all 12004) 



Table 4: The number of DRGs sample 



K,s bin center 


J 


- K s > 2.3 




N 


N corr Error 


(1) 


(2) 


(3) (4) 


17.35 


53 


53 19 


17.85 


120 


120 30 


18.35 


250 


270 50 



Note. — Description of the columns are the same as in the footnote of Table [T] 



Table 5: Adopted SED model and local LF parameters in the number count for early-type 
galaxies models 



Metalicity 


Star Fomation Rate 




0* a 


Z/Z 


Gyr 




Mpc~ 3 x 10~ 3 


(1) 


(2) 


(3) 


(4) (5) 


2.5 


exponential r = 0.3 


-24.0 


1.3 -0.85 



Note. - 



(2) Star formation and the e-folding time assuming Pure Luminosity Evolution (PLE). 



Table 6: Adopted local LF and SED parameters for star- forming galaxy number count models 



Model SED parameters 
Bolometric Luminosity starburst age IMF © Extinction Curve 



L/L 


Gyr 








(1) 


(2) 


(3) 


(4) 


(5) 


>10 10 


0.6 


1.35 


1.0 


SMC 


>10 n 


0.6 


1.35 


0.9 


SMC 


>10 12 


0.5 


1.35 


0.3 


SMC 



Model luminosity function parameters 

submm 

lg(L/L ) Mpc- 3 x 10~ 3 
Serjeant & Harrison (2005) # 7.89 2.5 0.458 3.36 



Note. - 

(1) L>1O 12 L SED is the model fit for the archetypical dERO HR10 

(3) Salpeter IMF is assumed 

(4) 6 is the opening angle of the star forming region 

(5) Small Magellanic Cloud extinction curve is assumed 

(#) Converted from H = 65 kms^Mpc -1 in Serjeant & Harrison (2005) to H = 70 kms^Mpc -1 



Table 7: The parameters for the evolutionary models for dEROs described in the text 



Population 


evolution to z=l 


evolution to z=3 




h 


9i 




92 


(1) 


(2) 


(3) 


(4) 


(5) 


STFG 


3.2 


0.9 


2.2 


0.4 


LIRG 


2.5 


3.2 


2.2 


2.7 


ULIRG 


2.7 


3.3 


2.2 


2.9 



